Molecular markers for the assessment of genetic variability in threatened plant species

La variabilità genetica è un importante strumento per lo studio e la conservazione della biodiversità in specie rare e minacciate di estinzione. Durante il mio dottorato mi sono quindi occupata di mettere a punto diverse metodologie molecolari al fine di valutare la diversità genetica in due specie rare della flora italiana che presentano problematiche diverse e specifiche. I marcatori arbitrari RAPD e i marcatori semi-arbitrari ISSR sono stati utilizzati per valutare la diversità genetica in Quercus crenata Lam. e per confermare l’ipotesi della sua origine ibridogena dalle due specie presunte parentali Quercus cerris L. e Quercus suber L., essendo Q. crenata presente in Italia settentrionale dove Q. suber è attualmente assente. I marcatori SSR o microsatelliti sono invece stati messi a punto su una specie a rischio di estinzione, endemica dell’Appennino Tosco-Emiliano, Primula apennina Widmer, applicando una metodologia specifica, basata sulla costruzione di una libreria genomica arricchita per l’isolamento di primer specifici. I marcatori RAPD e ISSR, utilizzati su un totale di 85 campioni, hanno mostrato alti livelli di diversità molecolare entro le specie studiate, eccetto per Q. suber le cui popolazioni rappresentano il margine orientale di distribuzione della specie, per questo più sottoposte ad impoverimento genetico. Oltre alla cluster analysis (UPGMA) e alla Analisi delle Componenti Principali effettuate per entrambi i marcatori, che confermano l’ipotesi dell’origine ibrida degli individui di Q. crenata diffusi in Italia Settentrionale, sono stati calcolati l’indice di ibridità basato sul maximum likelihood, che dimostra una introgressione asimmetrica di Q. crenata verso il parentale caratterizzato da superiorità demografica (Q. cerris) e il test di Mantel. Quest’ultimo ha permesso di confrontare i due marcatori RAPD e ISSR utilizzati ottenendo una bassa correlazione, a conferma del fatto che, amplificando tratti differenti del DNA nucleare, i dati non sono sovrapponibili, sebbene forniscano risultati analoghi. Per l’isolamento di loci microsatelliti ipervariabili ho utilizzato il protocolllo FIASCO (Fast isolation by AFLP of sequences containing repeats- Zane et al. 2002) che permette di costruire una libreria genomica arricchita partendo dal DNA estratto da P. apennina. Tale procedura ha previsto la digestione del DNA genomico per la produzione di una miscela di frammenti di DNA. Tramite ibridazione con opportune sonde sono stati isolati i frammenti contenenti i microsatelliti. Sequenziando i cloni ricombinanti, ho ottenuto sequenze contenenti repeats sulle cui regioni fiancheggianti sono stati costruiti 15 coppie di primer che potranno, in seguito, essere utilizzate per definire la quota di riproduzione clonale in P. apennina e per valutare la diversità genetica delle popolazioni che coprono l’areale di distribuzione della specie. Data la loro natura altamente variabile e la loro abbondanza nel DNA, gli SSR saranno, come i marcatori RAPD e gli ISSR, ugualmente validi per lo studio della variabilità genetica e per l’analisi di problematiche specifiche legate alle specie rare.

Genetica molecolare dell'autismo: studi di associazione ed analisi di geni candidati

Autism is a neurodevelpmental disorder characterized by impaired verbal communication, limited reciprocal social interaction, restricted interests and repetitive behaviours. Twin and family studies indicate a large genetic contribution to ASDs (Autism Spectrum Disorders). During my Ph.D. I have been involved in several projects in which I used different genetic approaches in order to identify susceptibility genes in autism on chromosomes 2, 7 and X: 1)High-density SNP association and CNV analysis of two Autism Susceptibility Loci. The International Molecular Genetic Study of Autism Consortium (IMGSAC) previously identified linkage loci on chromosomes 7 and 2, termed AUTS1 and AUTS5, respectively. In this study, we evaluated the patterns of linkage disequilibrium (LD) and the distribution of haplotype blocks, utilising data from the HapMap project, across the two strongest peaks of linkage on chromosome 2 and 7. More than 3000 SNPs have been selected in each locus in all known genes, as well as SNPs in non-genic highly conserved sequences. All markers have been genotyped to perform a high-density association analysis and to explore copy number variation within these regions. The study sample consisted of 127 and 126 multiplex families, showing linkage to the AUTS1 and AUTS5 regions, respectively, and 188 gender-matched controls. Association and CNV analysis implicated several new genes, including IMMP2L and DOCK4 on chromosome 7 and ZNF533 and NOSTRIN on the chromosome 2. Particularly, my contribution to this project focused on the characterization of the best candidate gene in each locus: On the AUTS5 locus I carried out a transcript study of ZNF533 in different human tissues to verify which isoforms and start exons were expressed. High transcript variability and a new exon, never described before, has been identified in this analysis. Furthermore, I selected 31 probands for the risk haplotype and performed a mutation screen of all known exons in order to identify novel coding variants associated to autism. On the AUTS1 locus a duplication was detected in one multiplex family that was transmitted from father to an affected son. This duplication interrupts two genes: IMMP2L and DOCK4 and warranted further analysis. Thus, I performed a screening of the cohort of IMGSAC collection (285 multiplex families), using a QMPSF assay (Quantitative Multiplex PCR of Short fluorescent Fragments) to analyse if CNVs in this genic region segregate with autism phenotype and compare their frequency with a sample of 475 UK controls. Evidence for a role of DOCK4 in autism susceptibility was supported by independent replication of association at rs2217262 and the finding of a deletion segregating in a sib-pair family. 2)Analysis of X chromosome inactivation. Skewed X chromosome inactivation (XCI) is observed in females carrying gene mutations involved in several X-linked syndromes. We aimed to estimate the role of X-linked genes in ASD susceptibility by ascertaining the XCI pattern in a sample of 543 informative mothers of children with ASD and in a sample of 164 affected girls. The study sample included families from different european consortia. I analysed the XCI inactivation pattern in a sample of italian mothers from singletons families with ASD and also a control groups (144 adult females and 40 young females). We observed no significant excess of skewed XCI in families with ASD. Interestingly, two mothers and one girl carrying known mutations in X-linked genes (NLGN3, ATRX, MECP2) showed highly skewed XCI, suggesting that ascertainment of XCI could reveal families with X-linked mutations. Linkage analysis was carried out in the subgroup of multiplex families with skewed XCI (≥80:20) and a modest increased allele sharing was obtained in the Xq27-Xq28 region, with a peak Z score of 1.75 close to rs719489. In this region FMR1 and MECP2 have been associated in some cases with austim and therefore represent candidates for the disorder. I performed a mutation screen of MECP2 in 33 unrelated probands from IMGSAC and italian families, showing XCI skewness. Recently, Xq28 duplications including MECP2, have been identified in families with MR, with asymptomatic carrier females showing extreme (>85%) skewing of XCI. For these reason I used the sample of probands from X-skewed families to perform CNV analysis by Real-time quantitative PCR. No duplications have been found in our sample. I have also confirmed all data using as alternative method the MLPA assay (Multiplex Ligation dependent Probe Amplification). 3)ASMT as functional candidate gene for autism. Recently, a possible involvement of the acetylserotonin O-methyltransferase (ASMT) gene in susceptibility to ASDs has been reported: mutation screening of the ASMT gene in 250 individuals from the PARIS collection revealed several rare variants with a likely functional role; Moreover, significant association was reported for two SNPs (rs4446909 and rs5989681) located in one of the two alternative promoters of the gene. To further investigate these findings, I carried out a replication study using a sample of 263 affected individuals from the IMGSAC collection and 390 control individuals. Several rare mutations were identified, including the splice site mutation IVS5+2T>C and the L326F substitution previously reported by Melke et al (2007), but the same rare variants have been found also in control individuals in our study. Interestingly, a new R319X stop mutation was found in a single autism proband of Italian origin and is absent from the entire control sample. Furthermore, no replication has been found in our case-control study typing the SNPs on the ASMT promoter B.

Harnessing the cellular and molecular complexity of high-risk neuroblastoma to investigate novel potential treatment approaches

Neuroblastoma (NB) is the most common type of tumor in infants and the third most common cancer in children. Current clinical practices employ a variety of strategies for NB treatment, ranging from standard chemotherapy to immunotherapy. Due to a lack of knowledge about the molecular mechanisms underlying the disease's onset, aggressive phenotype, and therapeutic resistance, these approaches are ineffective in the majority of instances. MYCN amplification is one of the most well-known genetic alterations associated with high risk in NB. The following work is divided into three sections and aims to provide new insights into the biology of NB and hypothetical new treatment strategies. First, we identified RUNX1T1 as a key gene involved in MYCN-driven NB onset in a transgenic mouse model. Our results suggested that that RUNX1T1 may recruit the Co-REST complex on target genes that regulate the differentiation of NB cells and that the interaction with RCOR3 is essential. Second, we provided insights into the role of MYCN in dysregulating the CDK/RB/E2F pathway controlling the G1/S transition of the cell cycle. We found that RB is dispensable in regulating MYCN amplified NB's cell cycle, providing the rationale for using cyclin/CDK complexes inhibitors in NBs carrying MYCN amplification and relatively high levels of RB1 expression. Third, we generated an M13 bacteriophage platform to target GD2-expressing cells in NB. Here, we generated a recombinant M13 phage capable of binding GD2-expressing cells selectively (M13GD2). Our results showed that M13GD2 chemically conjugated with the photosensitizer ECB04 preserves the retargeting capability, inducing cell death even at picomolar concentrations upon light irradiation. These results provided proof of concept for M13 phage employment in targeted photodynamic therapy for NB, an exciting strategy to overcome resistance to classical immunotherapy.

Utilizzo integrato di metodiche molecolari per la ricerca e la caratterizzazione genetica degli agenti infettivi in carnivori domestici e selvatici

Nel periodo compreso tra il 2019 e il 2022 sono state testate differenti matrici biologiche di carnivori domestici e selvatici provenienti dall’Italia e da altri Paesi europei (Norvegia, Romania). Diversi saggi molecolari, tra cui real-time PCR, end-point PCR, semi-nested PCR, retrotrascrizione e rolling circle amplification, sono stati utilizzati per ricercare il DNA o l’RNA genomico di virus e batteri. Il sequenziamento dell’intero genoma o di geni informativi dei patogeni identificati ne ha inoltre consentito la caratterizzazione genetica e l’analisi filogenetica. Gli studi, svolti presso il Dipartimento di Scienze Mediche Veterinarie dell’Università di Bologna, erano focalizzati nei confronti di alcuni virus a DNA, come Carnivore protoparvovirus 1 in lupi dall’appennino italiano e cani dalla Romania, adenovirus canino di tipo 1 e 2 in cani e lupi provenienti dal territorio nazionale, circovirus canino in cani e lupi italiani e volpi rosse e artiche della Norvegia; virus a RNA, come il canine distemper virus in faine recuperate nel territorio italiano e il calicivirus felino in gatti con diagnosi di poliartrite; e batteri appartenenti alla specie Anaplasma phagocytophilum in gatti deceduti e sottoposti a necroscopia in Italia. Dai risultati ottenuti è emerso che gli agenti infettivi indagati circolano nelle popolazioni di carnivori domestici e selvatici in forma asintomatica o determinando talvolta sintomatologia clinica. In alcuni animali testati è stata rilevata la coinfezione con diversi agenti patogeni, condizione che può predisporre ad un aggravamento della sintomatologia clinica. Dall’analisi filogenetica sono emerse relazioni tra gli agenti infettivi rilevati nelle differenti specie animali suggerendone la trasmissione tra ospiti domestici e selvatici e confermando il ruolo epidemiologico svolto dei carnivori selvatici nel mantenimento dei patogeni nel territorio. Alla luce dei dati ottenuti, è importante sottolineare l’importanza delle misure di profilassi, in particolare la vaccinazione degli animali da compagnia, per ridurre la trasmissione e la diffusione degli agenti infettivi.

Gwas analysis for the identification of molecular basis responsible for yield related traits and disease resistance in durum wheat

The PhD thesis was developed in the framework of Innovar H2020 project. This project aimed at using genomics, transcriptomics and phenotyping techniques to update varietal registration procedure used in Europe for Value of Cultivation and Use (VCU) and Distinctiness Uniformity and Stability (DUS) protocols. The phenotypic and genotypic diversity of a durum wheat panel were assessed for different agronomic traits, connected with wheat development, disease resistance and spike fertility. A panel of 253 durum wheat varieties was characterized for VCU and DUS traits and genotyped with Illumina 90K SNP Chip array (Wang et al., 2014). GWAS analysis was performed, detecting strong QTLs confirmed also by literature review. Candidate genes were identified for each trait and molecular markers will be developed to be used for marker assisted selection in breeding programs. As for disease resistance, the panel was evaluated for resistance to Soil-Borne-Cereal-Mosaic-Virus (SBCMV). A major QTL, sbm2, was detected on chromosome 2B responsible for durum wheat resistance (Maccaferri et al., 2011). The sbm2 interval was explored by fine mapping on segregant population using KASP markers and by RNASeq analysis, detecting candidate genes involved in plant-pathogen reaction. As regards yield related traits, detailed analysis was performed on the GNI-2A QTL (Milner et al., 2016), responsible for increased number spike fertility. Fine mapping analysis was performed on durum panel identifying hox2 a strong candidate gene, codifying for transcription factor protein. The gene is paralogue of GNI-1 (Sakuma et al., 2019), and it has a 4 kbp deletion responsible for increased number of florets per spikelet. To conclude, the herein reported thesis shows a complete characterization of agronomic and disease resistance traits in modern durum wheat varieties. The results obtained will augment available information for each variety, identifying informative molecular markers for breeding purposes and QTLs/candidate genes responsible for different agronomic traits.

Integrating structural variant calling, annotation and prioritization into whole genome analysis workflows: a practical application in the molecular diagnosis of neurodevelopmental disorders

Background: WGS is increasingly used as a first-line diagnostic test for patients with rare genetic diseases such as neurodevelopmental disorders (NDD). Clinical applications require a robust infrastructure to support processing, storage and analysis of WGS data. The identification and interpretation of SVs from WGS data also needs to be improved. Finally, there is a need for a prioritization system that enables downstream clinical analysis and facilitates data interpretation. Here, we present the results of a clinical application of WGS in a cohort of patients with NDD. Methods: We developed highly portable workflows for processing WGS data, including alignment, quality control, and variant calling of SNVs and SVs. A benchmark analysis of state-of-the-art SV detection tools was performed to select the most accurate combination for SV calling. A gene-based prioritization system was also implemented to support variant interpretation. Results: Using a benchmark analysis, we selected the most accurate combination of tools to improve SV detection from WGS data and build a dedicated pipeline. Our workflows were used to process WGS data from 77 NDD patient-parent families. The prioritization system supported downstream analysis and enabled molecular diagnosis in 32% of patients, 25% of which were SVs and suggested a potential diagnosis in 20% of patients, requiring further investigation to achieve diagnostic certainty. Conclusion: Our data suggest that the integration of SNVs and SVs is a main factor that increases diagnostic yield by WGS and show that the adoption of a dedicated pipeline improves the process of variant detection and interpretation.